Magnetic properties of FeSe superconductor 
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A detailed magnetization study for the novel FeSe superconductor is carried out to investigate 
the behavior of the intrinsic magnetic susceptibility \ in the normal state with temperature and 
under hydrostatic pressure. The temperature dependencies of x and its anisotropy Ax = X|] — X± 
are measured for FeSe single crystals in the temperature range 4.2 — 300 K, and a substantial 
growth of susceptibility with temperature is revealed. The observed anisotropy Ax is very large 
and comparable with the averaged susceptibility at low temperatures. For a polycrystalline sample 
of FeSe, a significant pressure effect on x is determined to be essentially dependent on temperature. 
Ah initio calculations of the pressure dependent electronic structure and magnetic susceptibility 
indicate that FeSe is close to magnetic instability with dominating enhanced spin paramagnetism. 
The calculated paramagnetic susceptibility exhibits a strong dependence on the unit cell volume and 
especially on the height Z of chalcogen species from the Fe plane. The change of Z under pressure 
determines a large positive pressure effect on x which is observed at low temperatures. It is shown 
that the literature experimental data on the strong and nonmonotonic pressure dependence of the 
superconducting transition temperature in FeSe correlate qualitatively with calculated behavior of 
the density of electronic states at the Fermi level. 

PACS numbers: 74.70.Xa, 74.62.Fj, 75.10.Lp, 75.30.Cr 



I. INTRODUCTION 

Soon after the discovery of superconductivity in 
LaFeAsOi-^Fy, the superconductivity was also detected 
in the binary compound FeSei_a~ with transition tem- 
perature T c ~ 8 K. This compound possesses the sim- 
plest crystal structure among the new families of Fe- 
based superconductors, and consists of a stack of Fc 
square planar layers, which are tetrahedrally coordinated 
by Se atoms. Also the large pressure effect on transition 
temperature was later observed 2-4 with T c w 37 K at 
pressures P « 9GPa, indicating that FeSei-z is actu- 
ally a high temperature superconductor. Therefore, the 
superconducting FeSei-jc compound has attracted con- 
siderable attention and is a subject of intensive studies 
for the last years^— . The structural simplicity of FeSe 
favors experimental and theoretical studies of chemical 
substitution and high pressure effects, which are aimed 
at promoting a better understanding of a mechanism of 
the superconductivity. 

Upon cooling below room temperature, the tetragonal 
PA/nmm phase of FeSe i -a; undergoes a subtle distortion 
to the lower symmetry orthorhombic Cmma phase 3 -^— . 
This transition occurs within a broad temperature range, 
about 70-^100 K, depending on the stoichiometry of the 
samples. It was also found that the tetragonal phase un- 
dergoes structural transitions under high pressures (P > 
10 GPa) to the hexagonal non-superconducting PQ^mmc 
NiAs-type phase^i 3 -^, and then to its orthorhombic mod- 
ification (Pbnm, MnP-type)^ ' 13 ' 14 . The recent theoret- 
ical examination of stability regions of the high pres- 



sure phases of FeSei- a, 15 ' 16 indicated a possibility of 
metallization and superconductivity in the orthorhombic 
Pbnm phase under high pressure. 

Though a substantial increase of T c was clearly ob- 
served under pressur e - 2 ' 3 ' 17 " — , these studies did not de- 
tect any magnetic ordering. However, recent NMR stud- 
ies provided some indication of an incipient magnetic 
phase transition under pressure^. Recently, a static 
magnetic ordering has been detected above P ~ 1 GPa 
by means of zero-field muon spin rotation (ZF ^tSR) 21 ' 22 
and neutron diffraction 2 ^. These studies have revealed 
that as soon as magnetic ordering emerges, the magnetic 
and superconducting states apparently coexist, and both 
the magnetic ordering temperature Tn and T c increase 
simultaneously with increasing pressure. Also, it was re- 
cently found, that upon applying pressure the increase 
of T c in FeSei-z appeared to be nonmonotonic and ex- 
hibits a local maximum at P ~ 0.8 GPa, which is followed 
by a local minimum at P ~ 1.2 GPa 1 -^^^ . There- 
fore, there is still a considerable controversy regarding an 
interplay between electronic structure, magnetism and 
superconductivity in the FeSei-^ compound, especially 
under pressure. 

The experimental data on magnetic susceptibility of 
FeSei-a; system in the normal state are incomplete and 
contradicting^ 1 ^. Also, these data are mostly obtained 
on polycrystalline samples and often distorted by the 
presence of secondary magnetic phases of iron. In or- 
der to shed more light on the relation between magnetic 
and superconducting properties, it is very important to 
elucidate the intrinsic susceptibility of FeSei-a supercon- 
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ductors and investigate its evolution with doping, tem- 
perature, and pressure. Here we report on results of the 
experimental studies of magnetic properties for single- 
crystalline and polycrystalline FeSe samples of high qual- 
ity in the normal state. These studies include mea- 
surements of the temperature dependence of magnetic 
susceptibility and its anisotropy as well as the hydro- 
static pressure effects. The experimental investigations 
are supplemented by ab initio calculations of the elec- 
tronic structure and magnetic susceptibility of FeSe. The 
calculations are based on the local-spin-density approxi- 
mation (LSDA) of the density- functional theory (DFT). 
The results of experiments and calculations are used to 
analyze the nature of magnetism in FeSe and the basic 
mechanisms of its strong pressure dependence. 



II. EXPERIMENTAL DETAILS AND RESULTS 

The plate-like single crystals of FcSci_j. superconduc- 
tor were grown in evacuated quartz ampoules using the 
KCI/AICI3 flux technique with a constant temperature 
gradient of 5° C per cm along the ampoule length (tem- 
perature of the hot end was kept at 427° C, temper- 
ature of the cold end was about 380° C). Typical di- 
mensions of the produced single crystalline samples are 
(2-j-3) x (2 -=-3) x (0.3H-0.5) mm 3 . The tetragonal PA/nmm 
structure was demonstrated at room temperature by 
X-ray diffraction technique. The energy dispersive X- 
ray spectroscopy, performed on a CAMECA SX100 (15 
keV) analytical scanning electron microscope, revealed 
Fe:Se=l:0.96±0.02 composition denoted in the following 
as FeSe for simplicity. 

The study of magnetic properties of FeSe samples at 
ambient pressure was carried out at T = 4.2 -=- 300 K 
by using a SQUID magnetometer. The superconducting 
transition was detected within 6-^8 K. The magnetiza- 
tion dependencies M(H) in magnetic fields up to 5 T 
appeared to be close to linear, indicating that the concen- 
trations of ferromagnetic impurities are negligibly small. 
A typical temperature dependence of the magnetic sus- 
ceptibility for single crystalline FeSe samples is shown in 
Fig. [U As is seen, a substantial growth of susceptibility 
with temperature was revealed in the normal state, as 
well as a large magnetic anisotropy. 

The study of magnetic susceptibility of FeSe under he- 
lium gas pressure P up to 0.2 GPa was performed at fixed 
temperatures 78 and 300 K by using a pendulum-type 
magnetometer placed directly in the nonmagnetic pres- 
sure cell22 In order to measure the pressure effect with a 
reasonable accuracy, a sufficiently large mass of the sam- 
ple is required. We have used the FeSe sample, further 
called as "polycrystalline" FeSe, which was prepared by 
compacting a number of about 50 of small arbitrarily ori- 
ented single crystals inside of an aluminum foil cylinder. 
The total mass of the sample was about 200 mg. The 
measurements were carried out in the field H = 1.7 T 
and their relative errors did not exceed 0.5%. The exper- 
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FIG. 1: Temperature dependence of the magnetic susceptibil- 
ity in the normal state for the single-crystalline FeSe measured 
in the field H — 50 mT. Data corresponding to the magnetic 
field directions H _L c axis and H || c are denoted by O and 
A symbols, respectively. 



imental pressure dependencies x(P) a t different temper- 
atures are shown in Fig. [5J which demonstrates their lin- 
ear character. The negative value of the pressure effect, 
dlnx/d-P ~ —6.5 x 10~ 2 GPa -1 , was observed at tem- 
perature 300 K, whereas at T = 78 K the effect appeared 
to be positive, dlnx/dP ~ 10 x 10 -2 GPa -1 . The avail- 
able experimental and theoretical results on dlnx/dP for 
FeSe are compiled in Table U together with corresponding 
data for the relative FeTe compound for comparison. 
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FIG. 2: Pressure dependencies of the magnetic susceptibility, 
normalized to its value at P = 0, for the polycrystalline FeSe 
compound at temperatures 78 and 300 K. The solid lines are 
guides for the eye. 
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TABLE I: Pressure derivatives of the magnetic susceptibil- 
ity dlnx/dP (in units 1(T 2 GPa -1 ) for FeSe and FeTe com- 
pounds at different temperatures. 
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a From NMR Knight shift data of Ref. YK 
b Results for FeTe are taken from Ref. [24 



FIG. 3: Band structure of FeSe around the Fermi level (at 
eV) marked by a horizontal line. 



III. COMPUTATIONAL DETAILS AND 
RESULTS 

In order to analyze the magnetic properties of FeSe 
compound in the normal state, the ab initio calculations 
of the electronic structure and paramagnetic susceptibil- 
ity were carried out. At ambient conditions FeSe com- 
pound possesses the tetragonal PbO-type crystal struc- 
ture (space group PA/nmm), which is composed by alter- 
nating triple-layer slabs. Each iron layer is sandwiched 
between two nearest-neighbor layers of Se, which form 
edge-shared tetrahedra around the iron sites. The posi- 
tions of selenium layers are fixed by the structural pa- 
rameter Z, which represents the relative height of Sc 
atoms above the iron plane. The structural parameters 

of FeSe were determined by means of X-ray and neutron 
diffractionM 1 L9 1 i4 1 25. 

The main purpose of the present ab initio calculations 
was to evaluate the paramagnetic response in an external 
magnetic field and to elucidate the nature and features 
of magnetism in the normal state of FeSe compound. In 
the context of this task, the dependencies of the mag- 
netic susceptibility on volume, lattice parameters and 
temperature were addressed. Ab initio calculations of the 
electronic structure of FeSe were performed by employ- 
ing a full-potential all-electron relativistic linear muffin- 
tin orbital method (FP-LMTO, code RSPtSSiSI). The 
exchange-correlation potential was treated within the lo- 
cal spin density approximation^ of the DFT. The cal- 
culated basic features of electronic structure of FeSe ap- 
peared to be in a qualitative agreement with results of 
previous DFT calculations^— . 

As is seen in Figs. [3]and[H the calculated band struc- 
ture and density of electronic states (DOS) of FeSe in- 
dicate the presence of the hybridized predominantly d- 
like Fe electronic states close to the Fermi level E-p. 
The chalcogen p-states are situated well below Ep and 
slightly hybridized with the d-states of iron. Also, one 
can see a van Hove singularity in Fig. @] at ~ 40 meV 



below E-p. It should be noted, that a proximity of the 
van Hove singularity to the Fermi level is considered as 
the key ingredient for superconductivity in iron-based 
superconductors^. 
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FIG. 4: Total density of electronic states of FeSe (solid line) 
and the partial contribution of the iron d-states (dashed line). 
The Fermi level position at eV is marked by a vertical line. 

To evaluate the paramagnetic susceptibility of FeSe, 
the FP-LMTO-LSDA calculations of field-induced spin 
and orbital (Van Vleck) magnetic moments were carried 
out within the approach described in Ref. HH- The rela- 
tivistic effects, including spin-orbit coupling, were incor- 
porated, and the effect of an external magnetic field H 
was taken into account self-consistently by means of the 
Zeeman term: 

Hz = MbH- (2s + 1). (1) 

Here /iB is the Bohr magneton, s and 1 are the spin and 
orbital angular momentum operators, respectively. The 
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FIG. 5: Calculated paramagnetic susceptibility of FeSe as a 
function of the unit cell volume (in A 3 ). Z is taken to be 0.26. 
The arrows indicate the theoretical (1) and experimental (2) 
equilibrium volume values. 



FIG. 7: Calculated temperature dependence of the paramag- 
netic susceptibility of FeSe. Z is taken to be 0.26, the unit 
cell volume is fixed between the theoretical and experimental 
values. The dashed line is a guide for the eye. 
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FIG. 6: Calculated paramagnetic susceptibility of FeSe as a 
function of internal parameter Z for the experimental unit 
cell volume. The arrow indicates the experimental value of 
Z. 



field induced spin and orbital magnetic moments provide 
estimations of the related contributions to the magnetic 
susceptibility, Xspin and Xorb- For the tetragonal crystal 
structure the components of these contributions, Xi\\ an d 
X<j_j ar e derived from the magnetic moments calculated 
in the external field of 10 T, which was applied parallel 
and perpendicular to the c axis, respectively. 

It is found that magnetic response to the external field 
is very sensitive to the unit cell volume, as well as to 
the structural parameter Z ', which represents the rela- 
tive height of chalcogen species from the Fe plane. The 
calculated dependencies of susceptibility of FeSe as func- 
tions of the volume and parameter Z are given in Figs. 
[5] and [51 respectively. 



Also, the thermal effects are taken into account in or- 
der to calculate the temperature dependence of param- 
agnetic susceptibility for FeSe compound. In this case 
the field induced spin and orbital magnetic moments 
were evaluated by corresponding integration with the en- 
ergy derivative of the Fermi-Dirac distribution function 
f(E, /x, T), and the temperature dependence of x was ac ~ 
tually determined by taking into account the finite width 
of -df/dE (see Ref. [13 for details). It should be noted 
that the energy derivative of the Fermi-Dirac distribution 
corresponds to a Dirac delta function at low tempera- 
tures, having a sharp positive peak at the Fermi energy 
E-p. This steep behavior of —df/dE resulted in some 
instability in the numerical calculations of \i which are 
seen in Fig. [7] 



IV. DISCUSSION 

The experimental superconducting transition temper- 
atures 6-^8 K, obtained for the studied FeSei-a; samples, 
agree with those reported in literatures^—. Above T c , a 
substantial growth of susceptibility with temperature is 
revealed in the normal state of FeSe up to 300 K (Fig. [I]), 
what indicates the itinerant nature of electronic states of 
Fe at the Fermi energy. 

The total susceptibility in the absence of spontaneous 
magnetic ordering can be represented as the sum: 

Xtot = Xspin + Xorb + Xdia + XL, (2) 

where these terms correspond to the Pauli spin suscepti- 
bility, a generalization of the Van Vleck orbital paramag- 
netism, the Langevin diamagnetism of closed ion shells, 
and a generalization of Landau conduction electrons dia- 
magnetism, respectively. Obviously, the Xdia term does 



■5 



not provide any anisotropy, and for FcSe the Langevin 
diamagnetism of closed ion shells can be estimated ac- 
cording to Ref. |35| as Xdia — —0.2 x 10~ 4 emu/mol. 

In order to analyse the experimental data on x f° r 
FeSe we used the calculated paramagnetic contributions 
to susceptibility, Xspin and Xorb- It has been shown 3 - 4 - 
that for paramagnetic metallic systems the Stoner ap- 
proach underestimates the spin susceptibility, whereas 
the LSDA field-induced calculations take into account 
non-uniform induced magnetization density in the unit 
cell and provide more adequate description. The first- 
principles calculations of the paramagnetic susceptibility 
for FeSe revealed that this system is in close proximity 
to the magnetic critical point. This can be seen from a 
steep rise of x{V) an d x{%) i n Figs. [S]and|Hl respectively, 
above the corresponding experimental values of V and Z . 
In fact, the calculated Stoner enhancement S ~ 7 clearly 
indicates that FeSe is close to the ferromagnetic Stoner 
instability. 

For the unit cell volume chosen between the theoreti- 
cal and experimental values, the dominant spin contribu- 
tion to magnetic susceptibility of FeSe is estimated to be 
Xspin — 2.4 x 10~ 4 emu/mol. The averaged orbital Xorb 
term amounts to ~ 0.4 x 10 emu/mol), being of about 
15% of the total paramagnetism. From comparison of 
the calculated paramagnetic susceptibility for the ground 
state in Fig.0 Xpara = Xspin+Xorb ^ 2.8x 10~ 4 emu/mol, 
with the experimental value Xcxp at T — > K in Fig. [T] 
Xcxp = (X|| + 2xj_)/3 ~ 1.5 x 10 emu/mol, it is clear 
that calculated value Xpara has to be substantially com- 
pensated by a diamagnetic contribution in order to con- 
form with the experimental data. One can estimate the 
expected diamagnetic contribution to magnetic suscepti- 
bility of FeSe to be about (xox P — Xpara) ~ —1-3 x 10 -4 
emu/mol. This diamagnetism is comparable in absolute 
value with the paramagnetic contribution, being much 
larger than the estimated above Langevin diamagnetism 
Xdia- Apparently, it can be ascribed to the xl term in 
Eq. ©. 

According to the experimental data in Fig. [TJ the ob- 
served anisotropy of susceptibility Ax is large in FeSe, 
and even comparable with the averaged susceptibility it- 
self at low temperatures. It appears to be much larger 
than the calculated anisotropy of orbital contribution to 
paramagnetic susceptibility of FeSe, Axorb = Xorb|| — 
Xorbj_ — —0.1 x 10 -4 emu/mol. Therefore, in order to 
explain the experimental Ax one can assume the presence 
of a substantial and presumably anisotropic diamagnetic 
contribution from conduction electrons. 

To calculate the Landau diamagnetic contribution xl 
is a rather difficult problem 36,37 . The free-electron Lan- 
dau approximation, which is often used for estimations, 
provides Xl that equals — | of the Pauli spin suscepti- 
bility. However, for many metallic systems the diamag- 
netism of conduction electrons xl can be many times 
larger than the free-electron Landau estimate Xlj an< ^ 
such anomalous and anisotropic diamagnetism is often 
determined by the presence of quasi-degenerated states 



with small effective masses at Ep (see Ref. and ref- 
erences therein). As can be seen in Fig. [3l in FeSe the 
quasi-degenerate states with small effective masses exist 
at Ef around the symmetry points T and Z, where the 
band degeneracies are lifted by the spin-orbital coupling. 
Such band structure features are of particular impor- 
tance in connection with a manifestation of the anoma- 
lously large and anisotropic xl , which was found to orig- 
inate from the similar degeneracy points 37 . It should 
be emphasized that rigorous theoretical analysis of xl is 
a rather cumbersome procedure, which goes beyond the 
aims of the present work. At this stage, we have iden- 
tified appropriate electronic states near Ep as possible 
sources of the large and anisotropic conduction electrons 
diamagnetism in FeSe. 

The theoretically evaluated temperature dependence of 
the paramagnetic susceptibility Xpara 

(T) in Fig. [3 which 
takes into account the finite width of the energy deriva- 
tive of the Fermi-Dirac distribution function, actually 
provides only a slight increase in x with temperature. 
Thus, the observed substantial growth of x(T) is rather 
puzzling at present. It is presumably related to a fine 
structure of DOS at Ef, but one should expect that FeSe 
system is driven far from the ground state at room tem- 
peratures. At this stage we should admit that increase of 
the unit cell volume for the tetragonal PA/nmm phase 
(AV/V ~ 1% with temperature rising up to 300 K, see 
Ref. [12) can provide about 10% growth of paramagnetic 
susceptibility, according to the x(V) dependence in Fig. 
[5] However, this volume expansion does not explain the 
experimental x(T) i n Fig- HI Also, a change of the struc- 
tural parameter Z with temperature can be substantial 
and of importance due to the strong x(Z) dependence in 
Fig. [6l but to date the influence of temperature on Z was 
not studied in a systematic way. Therefore, noticeable 
temperature effects on the lattice parameters, chalcogen 
atom position Z , and electronic structure itself should be 
taken into account in a rigorous quantitative analysis of 
X(T) in FeSe. 

The measured pressure effects on magnetic susceptibil- 
ity of FeSe are intriguing and require a detailed examina- 
tion. Firstly, as can be seen in Table HI there is a striking 
sign difference for the pressure effects on x at low and 
room temperatures. Also, the absolute value of this ef- 
fect is substantially larger than that observed in strongly 
enhanced itinerant paramagnets 3 - 4 -, and appeared to be 
comparable with such pressure effect on x reported for 
the related FeTe compound 2 — (see Table [TJ) . 

It should be noted that the present experimental data 
on x(T 7 P) for FeSe are in reasonable agreement with 
the results of Ref. 20 on temperature and pressure de- 
pendencies of the NMR Knight shift K of FeSe in the 
normal state. As can be seen in Fig. [8j the tempera- 
ture dependence K(T) at ambient pressure reflects the 
corresponding dependence of the magnetic susceptibil- 
ity in Fig. Q] Assuming the latter to be governed by 
the spin susceptibility Xspin (T), the only temperature 
dependent contribution in K{T) can be determined as 
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FIG. 8: Temperature dependencies of the NMR Knight shift 
K in FeSe measured at ambient pressure (O) and at P = 
1.4 GPa (dashed line). The data are taken from Ref. [2p| . 
The inset shows dependence of K on the averaged magnetic 
susceptibility for FeSe, \ — (X|| + 2xx)/3, from Fig. Q] 



K spin {T) = ax spi „(T) with a ~ 5 x 10 2 %(emu/mol)- 1 , 
resulted from the slope of K vs \ linear dependence in 
inset of Fig. [5] In addition, this contribution is expected 
to be also responsible for the pressure effect on K. Us- 
ing a rough approximation Xspin(P) ~ x(T) we obtain an 
estimate of dlnx/dP ~ d \nK sp i n /dP, which is listed in 
Table |H 

In order to clarify the behavior of x(P) in FeSe, we 
carried out the ab initio calculations of paramagnetic sus- 
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FIG. 9: Pressure dependence of the internal chalcogen struc- 
tural parameter Z for FeSe. The solid line indicates the results 
of calculations from Ref. |32j. Experimental data on parame- 
ter Z in FeSe for the tetragonal phase at T = 190 K (A, Ref. 
7), T = 295 K (□, Ref. HI), T = 300 K (O, Ref. i), and for 
the orthorhombic phase at T = 16 K (O, Ref. [J) and T = 50 
K (•, Ref. 0). In inset: data of Ref. on the expanded scale. 



ceptibility as a function of pressure. These calculations 
are based upon the pressure dependent structural param- 
eters, which have been calculated and listed in Ref. |32|. 
As is seen in Fig. [9j the calculated behavior of Z(P) 
follows the available experimental data^^ ' 14 i 25 . In the 
course of corresponding calculations of x(P) f° r FeSe we 
evaluated the pressure derivative dln^/dP ~ 8 x 1CP 2 
GPa -1 in the range 0-^1 GPa. The evaluated deriva- 
tive appeared to be in a qualitative agreement with the 
experimental low temperature data (see Table H]). 

As part of these electronic structure calculations for 
FeSe, we also obtained a pressure dependence of the den- 
sity of states at the Fermi level, which is presented in 
Fig. [TUl In addition to the calculated structural param- 
eters from Ref. [32|, we also employed the small upward 
shift AZ = +0.004 to start the optimized Z(P) depen- 
dence from the experimental value Z — 0.26 (see Z(P) 
behavior in Fig. [9]). The corresponding two sets of cal- 
culations demonstrate in Fig. [TTJ]a tolerable variation of 
N(E-p) behaviors in FeSe under pressure, depending on 
a small adjustment of Z(P = 0) between the theoretical 
and experimental values. 

With the aim to elucidate the main mechanism of the 
experimentally observed strong increase of the magnetic 
susceptibility of FeSe under pressure at low temperatures, 
we have also analyzed the pressure effect in terms of the 
corresponding changes of the volume and Z parameters 
by using the relation: 



d lnx 
dP 



d lnx d XxiV 

X 



d\nV 



dP 



d lux dZ 
dZ X dP 



(3) 



The required values of the partial volume and Z deriva- 
tives of x can he estimated from the results of ab ini- 
tio calculations presented in Figs. [5] and [51 and were 
found to be d\nx/dlnV ~ 8 and dlnx/dZ ~ 65 for 
the values of V and Z close to the experimental data. 
The optimized value dlnV/dP = -3 x 10" 2 GPa" 1 is 
taken for the compressibility of FeSe 32 . This calculated 
compressibility agrees closely with the experimental val- 
ues -3.1 ± 0.1 x 10~ 2 GPa" 1 iiii 3 .. Also, the optimized 
value dZ/dP ~ 0.55 x 10 -2 GPa -1 was adopted for eval- 
uation of Eq. ([3]). As is seen in Fig. [5J this value of 
dZ/dP is in agreement with the experimental data at 
low pressures. As far as all parameters entering Eq. ([3]) 
are estimated, the first term in ([3]) results in a large nega- 
tive value dlnx/dhiy x dlnV/dP ~ -24 x 10~ 2 GPa" 1 , 
whereas the second term appears to be large and pos- 
itive: dlnx/dZ x dZ/dP ~ 36 x 10~ 2 GPa" 1 . The 
both terms in Eq. © taken together yield the estimate 
dlnx/dP ~ 12 x 10~ 2 GPa -1 for FeSe, which is consis- 
tent with the low temperature experimental data. This 
estimate is also close to the ab initio calculated pressure 
derivative based on theoretical lattice parameters from 
Ref. HI (dlnx/dP ~ 8 x 10~ 2 GPa" 1 , see Table HJ). Actu- 
ally, the difference in these evaluated values of dlnx/dP, 
depending on whether experimental or theoretical lattice 
parameters are employed, covers a reasonable range for 
the expected pressure effect on x- 
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P (GPa) 

FIG. 10: Calculated pressure dependencies of the density of 
states at the Fermi level for FeSe.The pressure dependent 
structural parameters, including the lattice constants and 
chalcogen atom position Z, were taken from the calculations 
of Ref. 32 (O). For another set of N(E F ) calculations (□) 
the small upward shift AZ = +0.004 was employed to start 
from the experimental value Z — 0.26 (see Z(P) behavior in 
Fig- O • The solid lines are guides for the eye. 



Based on the results of calculations, presented in Figs. 
[5] and [6j the observed hydrostatic pressure effect on \ m 
FeSe at low temperatures can be represented as a sum 
of two large in size and competing contributions, related 
to the pressure dependence of the structural parameters 
V and Z . As a result, the experimental positive value of 
the pressure effect, dlnx/dP ~ 10 x 10 -2 GPa -1 , is de- 
termined by a dominating contribution from the change 
of Z under pressure. 

Actually, the nature of the large positive pressure ef- 
fect on x m FeSe is similar to that reported for FeTe 
compound 24 . However, in the case of FeTe such effect 
is twice more pronounced, and also takes place at room 
temperatures, whereas for FeSe dln^/dP is found to be 
negative at 300 K (see Table A} . The reason of this dif- 
ference is unclear and have to be elucidated. At the 
present stage one can presume, that the negative sign 
of dlnx/dP derivative is probably related to a nature 
of the observed anomalous growth of x(T) up to room 
temperatures (Fig. [JJ , which is not the case for FeTe 2 ^. 
It appears that at higher temperatures this anomalous 
growth of x(T) in FeSe is apparently reduced by applied 
pressure. 

Basically, the observed positive pressure effect on \ 
in FeSe at low temperatures correlates with the calcu- 
lated increase of the density of states at the Fermi level 
N(E-p) at low pressures (see Fig. IT0"|) . At higher pres- 
sures one can see nonmonotonic variation of N(E-p) in 
Fig. [10] which clearly exhibits consecutive maximum and 
minimum. It was recently shown^, that superconduct- 
ing transition temperatures T c of a number of iron-based 
superconductors correlate with the corresponding values 



of the density of states at the Fermi level, thus support- 
ing the BCS-like pairing mechanism in these systems. 
Remarkably, that the presently calculated behavior of 
N(Ep) under pressure (the upper curve in Fig. QUI with 
maximum at 0.5 GPa and minimum at 2.2 GPa) is qual- 
itatively consistent with the reported experimental de- 
pendencies of T C (P) in FeSe (corresponding maximum 
and minimum of T C (P) were observed at P ~ 0.8 GPa 
and P ~ 1.2 GPa, respectively^^ ) . The calculated 
pressure dependence of DOS at the Fermi level for FeSe 
with the structural parameters taken from Ref. [32| (the 
lower curve in Fig. I10|) also contains consecutive max- 
imum and minimum of N(Ep), which are substantially 
shifted to higher pressures. 



V. CONCLUSIONS 

The magnetic susceptibility of FeSe compound is found 
to rise substantially with temperature, which appar- 
ently points to the itinerant nature of the electronic 
states of Fe. The calculated paramagnetic susceptibil- 
ity Xpara(P) describes qualitatively the experimental de- 
pendence x(T), however the origin of the observed about 
twofold increase of x U P to 300 K is puzzling. From 
comparison of the experimental values of susceptibility 
and its anisotropy with that calculated for paramagnetic 
contributions to x m FeSe, the additional anisotropic dia- 
magnetism is expected to be of the order of —1 x 10~ 4 
emu/mol, which can relate to the diamagnetism of con- 
duction electrons and presumably has its origin in the 
quasi-degenerate electronic states close to Ep. 

The measurements of magnetic susceptibility under 
hydrostatic pressure revealed a strong positive effect at 
low temperatures. This effect appeared to be compara- 
ble with that reported for the related FeTe compound 24 , 
whereas at room temperature the pressure effect for FeSe 
is found to be also strong, but negative. 

Our calculations indicate that paramagnetic suscepti- 
bility of the FeSe compound is substantially dependent 
on the unit cell volume V and the relative height Z of Se 
species above the Fe plane. It is shown that the observed 
at low temperatures large positive pressure effect on \ 
is related to the strong sensitivity of the paramagnetic 
susceptibility to the parameter Z, which determines the 
dominant positive contribution. The grounds of the neg- 
ative sign of dhrx/dP derivative in FeSe at 300 K are 
unclear and probably linked to a nature of the observed 
anomalous growth of x(^)- At present one can state that 
at higher temperatures this anomalous growth of x(^) i n 
FeSe is apparently reduced by applied pressure. 

The present ab initio calculations have demonstrated 
that for FeSe compound the nonmonotonic behavior of 
superconducting transition temperature with pressure 
qualitatively correlates with the density of electronic 
states at the Fermi level. This indicates a possible real- 
ization of the BCS-like pairing mechanism in this system. 
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